Superfluidity of Bose-Einstein condensates in toroidal traps with nonlinear lattices 
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Superfluid properties of Bose-Einstein condensates (BEC) in toroidal quasi-one-dimensional traps 
are investigated in the presence of periodic scattering length modulations along the ring. The 
existence of several types of stable periodic waves, ranging from almost uniform to very fragmented 
chains of weakly interacting and equally spaced solitons, is demonstrated. We show that these 
waves may support persistent atomic currents and sound waves with spectra of Bogoliubov type. 
Fragmented condensates can be viewed as arrays of Josephson junctions and the current as a BEC 
manifestation of the dc- Josephson effect. The influence of linear defects on BEC superfluidity has 
been also investigated. We found that for subcritical velocities, linear defects that are static with 
respect to the lattice (while the condensate moves in respect to both the optical lattice and the 
defect) preserve the BEC superfluidity. 

PACS numbers: 00.00.AA 



I. INTRODUCTION 

The possibility of particles flows in multiply connected 
superfluid and superconducting systems which can last 
for extraordinary long, in principle infinitely long, times 
is certainly one of the most striking manifestation of 
quantum mechanics at the macroscopic level This 
fact, well known for conventional superconductors and 
superfluids, has been recently observed also in atomic 
Bose-Einstein condensates (BEC) confined in toroidal 
traps with persisting circulating currents. In ad- 

dition to the fundamental interest of this result as a di- 
rect manifestation of the BEC superfluidity [4|, matter 
waves in multiply connected geometries also represent 
ideal systems for exploring nonlinear properties in the 
presence of nonzero density backgrounds. In this con- 
text, solitary matter waves in one-dimensional (ID) and 
2D toroidal traps have been investigated and the exis- 
tence of a bandgap structure and of gap solitons have 
been reported in Ref. [5[. Possibility of management of 
solitons, including acceleration and localization of matter 
waves, was explored in Ref. Q . The density distributions 
in a ring as a function of the relation between the healing 
length and the trap length was investigated in Ref. [7|. 
Solitary waves have also been studied for quasi ID circu- 
lar troughs of radially periodic 2D potentials |8[- 

All these studies refer to the case of uniform or peri- 
odic linear potentials present along the circumference of 
the trap. Moreover, several experimental studies of BEC 
in linear optical lattices (OL) have shown the occurrence 
of interesting phenomena related to phase relationships 
between BEC wavefunctions in different sites, these in- 
cluding the existence of Josephson currents @ , the prolif- 
eration of vortices in the Berezinskii-Kosterlitz-Thouless 
regime of BECs in 2D OLs [l(| , vortex nucleation in ro- 



tating lattices of BECs fTTJ] , and so on. 

On the other hand, it is known that interesting phe- 
nomena also arise when a linear optical lattice is com- 
bined with a periodic spatial modulation of the nonlin- 
earity, also known as a nonlinear OL. This is particularly 
true for ID settings for which it has been shown that 
periodic spatial modulations of the scattering length, 
achieved with the help of the Feshbach resonance tech- 
nique 12], can induce long lived Bloch oscillations IR 



Rabi oscillations [14|, and dynamical localization 



gap solitons in the presence of accelerated linear periodic 
potentials. Periodic modulations of the nonlinearity have 
also been shown to be effective to generate solitons in ID 
OLs [l6[ and in random potentials [l7j to stabilize mul- 
tidimensional solitons [18| and to generate vortex rings 
in a highly controllable manner |19f (for fresh review of 
this very active field of investigation see Ref. 20] ) . 

We remark that toroidal traps are routinely created 
in laboratories with the aid of magnetic fields [2l[ and 
nonlinear OLs have also been recently experimentally re- 
alized for BEC in linear settings [13]. Methods to create 
trapping potentials for BECs by means of rapidly moving 
laser beams for a variety of geometries, including toroids, 
ring lattices and square lattices, have also been recently 
developed [23| . 

To our knowledge, nonlinear OLs with a multiply con- 
nected geometry have not yet been realized. Considering 
the experimental rapid progresses in the field, however, it 
is reasonable to expect them to be soon available. A non- 
linear toroidal lattice could indeed be created by optically 
induced Feshbach resonances [2jJ| using an all-optical trap 
with an horizontal sheet beam and a ring shaped vertical 
higher order Gauss-Laguerre beam [25[ whose intensity 
is periodically modulated along the ring due to phase in- 
terference (see, e.g., Ref. HH). This setting would be 
especially effective in a ID limit (for details of respective 
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reduction see, e.g., Ref.@) and would permit to study 
the combined effects of nonlinearity and periodicity on 
the BEC superfluidity, a problem scarcely investigated 
in the literature. 

The aim of this paper is to study superfluid properties 
of BECs in toroidal traps with strong radial confinement 
and periodic modulations of the scattering length (non- 
linearity) in the azimuthal direction in the framework of 
the mean field Gross-Pitaevskii equation (GPE). In par- 
ticular, we demonstrate the existence of several types of 
stable periodic waves, ranging from almost uniform to 
very fragmented chains of weakly interacting and equally 
spaced solitons, which can support persistent currents 
and permit sound waves with interesting physical prop- 
erties. The backgrounds that may exist in such traps can 
be classified in terms of bifurcation diagrams in the pa- 
rameter space {p, p) , where p is the chemical potential 
and p is the density of the condensate. We analyze the 
bifurcations patterns and the stability properties of peri- 
odic matter waves both with zero (stationary states) and 
nonzero particle currents. The spectral properties of the 
linear excitations (sound waves) propagating in periodic 
condensates are also investigated and the Bogoliubov 
type features of the spectrum demonstrated. In the pres- 
ence of a pronounced fragmentation of the condensate 
(i.e., in the soliton chain limit) we find that persistent 
currents of particles are expressed in the form p m i n sin#, 
where p m i„ is the minimum density of the condensate 
and is the phase difference of the order parameter of 
neighboring solitons. Quite interestingly, this expression 
is similar to the supercurrent expression of a Josephson 
junction in the zero voltage state (dc- Josephson effect), 
this suggesting the interpretation of fragmented conden- 
sates as arrays of Josephson junctions and the persistent 
current as a BEC analog of the dc- Josephson effect [13] ■ 
Superfluidity properties are also investigated in the pres- 
ence of localized linear defects perturbing the nonlinear 
OL and acting on the condensate as obstacle. Direct 
numerical simulations show that, in conformity with the 
Landau criterion and in agreement with our sound wave 
analysis, the BEC superfluidity is preserved in presence 
of linear defects which are static with respect to the non- 
linear OL and that move in respect to the condensate 
with subcritical velocities. 

The paper is organized as follows. In Sec.|TT]we formu- 
late the mathematical model describing the condensate 
and analyze the ground states properties. Families of 
periodic solitons in parameter space and their stability 
properties both for zero (stationary) and non zero car- 
rying current states are discussed in in Sees. II A and 
II B, respectively. The spectra of sound waves in the 
condensate and its Bogoliubov structure are numerically 
calculated and discussed in Sec. III. In Sec. IV the su- 
perfluid stability properties of BEC are investigated by 
direct numerical integrations of the GPE in the presence 
of a linear defect. Finally, in the Sec.V, the main results 
of the paper are summarized. 



II. MODEL EQUATION 

We consider the following dimensionless GP equation 

iipt = -i> xx + U(x)\i>\ 2 ilj + V(x)ip, (1) 

where U (x) is a periodical nonlinear potential that ac- 
counts for interatomic interactions and V(x) is a local- 
ized linear potential modeling an obstacle perturbing the 
flow of the ground state. Without loss of generality, the 
period of the nonlinear lattice is chosen to be tt, i.e. 
U(x + tt) — U(x) and the circumference of the trap is 
fixed to L = Mir, with M being a positive integer. The 
annular geometry obviously implies cyclic boundary con- 
ditions 

1>(x,t) = il>(x + L,t), (2) 

and permits the existence of current carrying states, with 
the current, J, defined as 

j = ~ (^*v* - w>t) (3) 

(here the asterisk denotes complex conjugation). The 
properties of these states are investigated in the following 
two subsections. 

A. Ground states with zero current 

We start from the the analysis of the stationary states 
having zero current: J = 0, assuming that no linear de- 
fect is added to the system: V(x) = 0. In this case 
the order parameter has a constant phase, and can be 
searched in the form ip = ty(x)e~ ttJ,t , where p is the chem- 
ical potential and is real and solves the stationary GP 
equation 

f i* = -* xx + U(x)* 3 , (4) 

We are particularly interested in nonlinear potentials 
changing the sign of the interaction, i.e. in the situa- 
tions in which the repulsive and attractive interactions 
alternate along the length L of the trap. We remark that 
for a priori fixed particular solutions ty(x) it is possible 
to construct exact potentials U(x) for Eq. ((H) using the 
" inverse engineering" technique explained in Ref . [28j . An 
example of such a potential is given by 

nq 2 snlqx, k)(3k sn(qx, k) — 2) 4K(k) 

u ( x ) = —a T~i / s . 1 = 

4 1 + Ksn(g:r, K) tt 

with sn denoting the standard Jacobi elliptic function 
and K(k) the elliptic integral. Notice that U(x) is 
parametrized by the elliptic modulus k £ [0, 1] and de- 
pending on the value of k, it represents either a sequence 
of nonlinear single wells (for < k < \/5/3 — 1) or a 
series of nonlinear double wells (for \/5/3 — 1 < k < 1). 
In this paper we shall restrict our study to the case 
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of single-well potentials, only. We also remark that in 
spite of the apparent complexity of the potential ([5]), it 
can be approximated very well by a Fourier series with 
only few harmonics. Indeed, even for the limiting value 
k = — 1 we obtain for the coefficients of the Fourier 

series U{x) = (U) + Y, k ^ c k e 2lkx , where (U) « 0.239, 
|c ±1 | w 0.36, |c ±2 | « 0.12, |c ±3 | w 0.016, |c ±4 | « 0.001, 
and so on. This implies that the suggested periodic mod- 
ulation of the nonlinearity can be simulated with accu- 
racy of about 5% by employing only two laser beams with 
the wavelengths n and 7r/2 in dimensionless units. 

An exact analytic solution of Eq. Q with V = and 
with the potential U(x) taken as in Eq.((5j), can be written 
as 



(a), 



* a = \Jl + KSii(qx, k), 



(6) 



with the corresponding chemical potential given by p a — 
n 2 q 2 /A. This solution is just a member of a family of pe- 
riodic solutions which are parameterized by the chemical 
potential p and identified by the period Lq. Note that 
although Lo does not need to coincide with the period of 
the nonlinear potential U(x), a relation L = mLg with to 
integer, must be obviously satisfied. Since we have fixed 
L = Mtt the set of possible periods of the solutions is 
given by all to, which are integer divisors of M, i.e., M/m 
is an integer and, respectively, Lq — (M/m)n = jm. 
In Fig. [T] we illustrate the bifurcation diagrams of these 
family solutions obtained numerically in the parameter 
space (p, p), where p = (\^f\ 2 ) is the linear density of the 
condensate |31| denoting the average of any L -periodic 
function f(x) as (/) = f(x)dx, we have found 

that the potential ([5]) satisfies (U) > 0, i.e. the poten- 
tial is repulsive in average. In the present paper we re- 
strict mainly to branches bifurcating from linear solutions 
for which the limiting transition p — > is possible. In 
this limit the ground state bifurcates from the uniform 
solution ty c = y/Ji/y/ (U) + 0(p 3 ^ 2 ) and the respective 
branch has the period 7r. The bifurcation points of the 
other solutions having linear limit can be designated by 
p m = Am 2 I'M 2 . In each of such points, however there 
are two emergent modes: cos- and sin-like modes: 



M - Mr, 



ir U{^)^{2x)dx 

M - Mm 

2 j;u(M±)^{2x)dx 



sin 



2to 
2m 



(7) 
(8) 



Note that, formally, the ground state can be obtained 
from Eq.([7]) by putting to = 0. 

From Figs. QJa) and Gib) we observe that an increase 
of the ground state density leads to the growth of the 
chemical potential, as is expected for a BEC with repul- 
sive interatomic interactions. However, when the number 
of particles increases the nonlinear interactions become 
important and the condensate becomes denser in the ar- 
eas of attractive interactions. At the point where the 
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FIG. 1: (Color online) [(a) and (b)] Families of static solutions 
for the potential ([5]) centered at x = tv/2 with k — 0.25 [solid 
black lines in panels (c)-(f)] bifurcating from linear modes for 
the cases: 

Lq — tt [red line (1)], Lq — 2n [black (2) and cyan (3) lines], 
and Lq = 4tt [purple (4) and green (5) lines], (c) The order 
parameter of the solutions is shown in panels (a) and (b) by 
red line (1). The solid red curve corresponds to the solution 
with p — —8 and the solution at the folding point p = 0.08 
is shown by the dashed blue line, (d) The order parameter 
of the solutions is shown in (a) and (b) by the black line (2) 
for p — — 8 (solid red line) and for p = 0.5 (dashed blue 
line), (e) The same as in panel (c) but for the solutions 
corresponding to p = 2.5 (dashed blue line) and p = —8 
(solid red line) of the family shown by cyan line (3). (f) The 

same as in panel (c) but for the solution corresponding to 
p — —8 (dashed blue line) of the family shown by the purple 
line (4) and for the solution corresponding to p = —8 (solid 
red line) of the family shown by the green line (5). 



bifurcation diagram turns back, referred below as fold- 
ing point, the chemical potential decreases [see the red 
line (1) in Figs. [TJa) andQJb)], and, for some average 
density, the chemical potential becomes negative. This 
happens because for larger densities most of the parti- 
cles concentrate in the areas with attractive interactions 
and, therefore, the whole condensate behaves as having a 
negative scattering length. Mathematically, this is clear 
from the relation 

pp = H = ^ 2 x ) + {U* i ), (9) 

from which we see that a negative p implies (C/^ 4 ) < 
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FIG. 2: (a) Solid lines show the dependencies of the lin- 
ear density on the chemical potential for the solutions with 
I/O < 4-7T and having no more than one bright soliton in each 
area of attractive interaction. The thick black and red lines 
correspond to the modes bifurcating from the linear solutions 
with the periods n and 2ir and, at large negative /i, having one 
soliton in each area of the attractive interactions. The other 
two modes originating from linear ones are shown by the green 
and purple lines. Thin blue lines show the dependencies for 
the solutions that do not have a linear limit. The dashed 
thick cyan line shows a bifurcation curve which corresponds 
to the solutions originating from a linear one and having two 
solitons in each area of the attractive interactions, (b) The 
field distributions for the solution corresponding to the points 
of the bifurcation curve denoted by bl and b2 in panel(a) 
(fi — —0.07) (black solid line). The red dashed lines show the 
solution with p, = —2.5 and belonging to the same bifurcation 
branches, (c) The same as in (b) but for the points marked 
in panel (a) by the circles cl and c2; the chemical potentials 
in this case are p = 0.125. 



(note, that in the folding point /i*, where dp,/ dp = we 
have the relations /i» = — = ^—)- Thus a condensate 
with sufficiently large density can be considered as an ar- 
ray of in-phase bright matter solitons located at the areas 
with negative scattering length; see Fig.QJc). Obviously, 
the interaction between the solitons decays exponentially 
as the density of the condensate increases, transforming 
in an array of in-phase matter solitons. 

On the other hand, from the diagram bifurcating from 
the point p = 1 [black (2) and blue (3) lines in Figs. Ufa) 
and [TJb)] we observe that the solution bifurcating from 
the sin mode [Eq.0] at large negative p can be seen as 
a chain of out-of-phase (...0 — 7r — — 7r...) bright solitons. 
The solutions belonging to the second branch bifurcating 
from Eq. ([8} in the strongly nonlinear limit also transform 
into a chain of pairs of bright solitons; in this case there 
are two matter solitons in each domain of attractive in- 
teractions. 



It is worth pointing out that the solution families 
shown in Fig. [T] reveal a peculiar feature, namely the 
kinetic energy of the condensate at the folding point is 
much larger than the potential energy ("F^) (U^ 4 ). 
We also note that not all periodical solutions bifurcate 
from linear modes: There are branches having no zero 
nonlinearity limit which is however, expectable in view 
of the analogous behavior of the localized modes [29| . 
Their bifurcation diagrams have folding points or they 
merge with other solutions. In Fig. [2] the bifurcation di- 
agrams are shown for all possible solutions which in the 
strongly nonlinear limit have no more than one soliton 
per period of the nonlinear potential and the period of 
these solutions is smaller than four periods of the nonlin- 
ear potential. 



B. Periodic matter waves with non zero current 

We now consider periodic states carrying nonzero cur- 
rents J ^ 0. Such states can be represented in the form 



*(a?) - ^fpjx)e l6{x \ 6{x) = / v(x)dx, (10) 



where v(x) is the superfluid velocity and p(x) the con- 
densate density. It is straightforward to show that such 
solutions exist only if J = v(x)p(x) = const [see the def- 
inition, Eq.Q], i.e. if dJ/dx — 0, with the density p 
satisfying the equation 

ld 2 p 1 fdpV J 2 TTf . 2 n , , 

(note that zero-current states are characterized by den- 
sities which are strictly positive so the above equation is 
nonsingular). In the following we consider Lo-periodic 
density distributions 



p(x + L ) = p(x) 



(12) 



of the form p(x) = p + pi(x), with p = (p{x)) the mean 
density and pi(x) a periodic variation of the density of 
period Lq. The condition J = const then implies that 
also the hydrodynamic velocity is of the form v(x) = vq + 
vi(x) with vq = (v) the average superfluid velocity and 
Vi a periodic function of period Lq. The order parameter 
^(x) has the form of a nonlinear Bloch state 



fv {%) 



(13) 



with the mean velocity vo playing the role of quasimo- 
mentum and with f Vo (x) being a complex periodic func- 
tion of period Lq. For the single valuedness or the wave- 
function ^(x), the average superfluid velocity must be 
quantized with respect to the length L of the ring ac- 
cording to 
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FIG. 3: (a) The current vs the average velocity for fi = —2 
(blue curve) and fj, = — 8 (red curve), (b) The average density 
vs the average velocity for fi = — 2. 



where p is an integer. We should remark that the analogy 
exploited above between superfluids in a ring and Bloch 
electrons in crystals was first discussed by Bloch [3(| and 
further developed in Ref.Q for BEC in rotating traps. 
Obviously, in the linear limit f vo (x) is simply reduced to 
y/~p, and, taking into account, that in this limit p — » 0, 
one can neglect the spatial variations of the density and 
conclude that the bifurcation curves start at p, = Vq. 

It is of interest to investigate how the condensate with 
a nonzero current behaves when its average velocity Vq 
is changed. This is done in Fig. |3] where the current 
and the average density versus the average velocity are 
depicted for two fixed values of the chemical potential. 
Note that the current depends on the average superfluid 
velocity periodically, acquiring maxima at semi-integer 
values Vo = p/2, while at integer values of Vo the cur- 
rent vanishes. This is a direct consequence of the Bloch 
form (fl~3|) of the wavefunction implying that nondegen- 
erate states with given \i and vq are equivalent to the 
state with the same chemical potential p but with the 
velocity vq + p, where p is an integer. Also note that at 
v o = the solution in Fig. [3] is carrying zero current [as 
the symmetric state belonging to the bifurcation branch 
(1)] while at vo = 1 the obtained periodic state is the 
asymmetric mode belonging to the branch (2) shown in 
Fig. [1] which is also a zero-current state. 

Another interesting property to remark is the quasi si- 
nusoidal dependence of the current J on vq for BEC wave- 
functions fragmented as chains of solitons on a nonzero 
background. The particle density of such condensates 
is depleted (enhanced) in correspondence with the re- 
gions of repulsive (attractive) interaction. Direct numer- 
ical simulations show that the stronger is the depletion 
of the condensate in the repulsive regions the better the 
current is approximated by a sinusoidal law. In the limit 
of very strong depletion (i.e., for a chain of strongly lo- 
calized and weakly overlapping solitons), it becomes 

J ~ ^\/ 3 MPmmSin6), (14) 

where p m in — — 8/V max((7) •sech 2 (- v /— Jj/k/2) denotes the 
minimum density of the condensate and 9 is the phase 
difference of the order parameter between neighboring 
solitons. We also note that in the case of strong localiza- 
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FIG. 4: (Color online) Temporal evolution of the BEC den- 
sity for a current close to the critical value. Red color corre- 
sponds to the maximum density and blue color corresponds 
to the minimum density. The width of the window adjusted 
to show four periods of the initial field distribution. The pa- 
rameters are jj, — —2, vo = -p=, and L = 6iir. 



tion the maximum condensate density corresponding to 
the soliton peaks, located at x = n/2 + 7rp where p is an 
integer, is given by p max = 2pJ min(U). 

Since the phase difference 9 is constant for adjacent 
solitons, the above expression of J is similar to the 
supercurrent flowing at the zero voltage in a conven- 
tional superconductive Josephson junction (e.g., to the 
dc-Josephson effect) (2?J ■ This suggests the interpreta- 
tion of the fragmented condensates as an array of Joseph- 
son junctions and the persistent current as a BEC mani- 
festation of the dc-Josephson effect [13] . The depleted ar- 
eas between solitons (i.e., where the interaction is repul- 
sive) play the role of insulating barriers of conventional 
Josephson junctions. Following this analogy we can say 
that for stable symmetric ground states the depleted ar- 
eas are the analogous of 0-junctions, while for the stable 
asymmetric states the depleted areas correspond to n- 
j unctions. Also note that the role of the electric potential 
across the junction is played by the difference of chemical 
potentials between neighboring solitons. Since the frag- 
mented state is a stationary state, it has a fixed chemical 
potential and, therefore, the flow of particles across the 
depleted regions (supercurrent) in the analogy also occurs 
in to the presence of zero voltage. Following this analogy, 
one can also expect that the nonzero-current stationary 
solutions become unstable and disappear when the cur- 
rent approaches a critical value (critical current). This is 
shown in Fig.2]where the temporal evolution of the BEC 
density is depicted for a current very close to its critical 
calue. We see that in this case the condensate becomes 
very sensitive to perturbations and develops instabilities 
that eventually lead to its destruction (see the next sec- 
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tion). We also observe that the critical current decays 
exponentially with the height and width of the effective 
nonlinear barrier, although it never vanishes completely. 



III. SOUND WAVES OF PERIODIC SOLITON 

The time evolution of sound waves against nonlinear 
periodic backgrounds of the form (fT3|) is governed by the 
linear equation 

iipt = -tp xx - 2iv (p x + (2U{x)\f Vo \ 2 + vl - fi)cp 

+u(x)f v y. (is) 

Looking for a monochromatic excitation, ip — 
A(x)e-^ u]+ ^ t + B* (x)e 4(w -^)* , where u(k) and j(k) are 
real, and splitting exponents of the respective ODEs one 
obtains the system of coupled equations with periodic 
coefficients: 

(ui + ij)A = -A xx - 2iv A x + 
(2U(x)\f V0 \ 2 + vl - fj)A + U(x)flB, (16) 

-(w + i~f)B = -B xx - 2iv B x + 
(2U(x)\f V0 1 2 + vl - fiB + U(x)f: 2 A. (17) 

These equations allow (due to the Floquet theorem) 
for solutions of the form ~ ipk(x)e lkx with (fk(x) a pe- 
riodic solution of Eq. (fT5]) [it is worth noting, that the 
period of the function f Vo is not necessarily equal to the 
period of the nonlinear potential U(x)] because the coef- 
ficients in Eas.fTS]) are 2U(x)\tp \ 2 and U(x)ipl, with the 
background state ipo a periodic function of x. [In this pa- 
per we restrict our consideration to periodic backgrounds 
only] . 

The eigenvalues and the corresponding eigenfunctions 
can be found numerically: uj(k) yields the phase and the 
group velocities of the sound waves, while j(k) defines 
the stability of the ground state, the instability corre- 
sponding to positive j(k) (note that since the problem 
is Hamiltonian the real part 7 can appear only in pairs 
± 7 ). 

Now let us turn to the stability of the backgrounds 
with zero current whose spectrum is symmetric with re- 
spect to the transformation k — > —k. First, we consider 
the symmetric modes corresponding to the curve 1 in 
Fig.QJb). In this case the period of the solution coincides 
with the period of the nonlinear potential and the width 
of the Brillouin zone of the linearized problem is equal 
to 2. The solutions belonging to the lower branch of the 
bifurcation diagram are stable. The typical spectrum of 
the linear excitation is shown in Fig.[SJa). The spectrum 
is pure real and is of Bogoliubov type so v g = d^Lo ^ at 
k = 0. At the folding point the background becomes un- 
stable, a fact that is confirmed by our spectral analysis. 
In particular, we find that the instability is generated by 
the merging of the lowest branches of the dispersion curve 
[i.e., the ones corresponding to lower values of |w(fc)|)], 
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FIG. 5: (Color online) Sound wave spectra for the symmetric 
mode corresponding to curve 1 in Fig. \V[b) at /1 = 0.04 (a) 
and at fi = —2 (b) on the lower and upper branches of the 
bifurcation curve, respectively. Panels (c) and (d) refer to 
sound wave spectra for backgrounds modes taken on curve 
2 of Fig. Ob) at chemical potentials /1 = and /j, = —2, 
correspondingly. The black and blue curves correspond to 
ui(k) (left vertical axis) and the red curves correspond to ~f(k) 
(the right vertical axis). 



as illustrated in Fig. [U[b) for the case /x = —2 for which 
branches are completely merged. 

Next we consider the antisymmetric backgrounds cor- 
responding to curve 2 in Fig.[IJb). In this case the coeffi- 
cients in the linearized Eq. (|15p have the period n equal to 
the period of the potential only if the state is described 
by a pure real function VP, i.e., if the background has zero 
current. In the general case the period is equal to 27T and 
the width of the Brillouin zone is equal to 1. The spectra 
of the linear excitations on antisymmetric backgrounds 
are shown in Figs.[5][c) and^d). 

The spectra shown in Fig. [5] reveal Bogoliubov struc- 
ture; however, they behave differently as far as the sta- 
bility concerns. Specifically, the states with relatively 
low particle density are unstable; see Figj5{c). It is seen 
that there are two pairs of merging modes, and the in- 
stability appears exactly at the wavevectors where two 
modes merge. On the other hand, at some threshold 
density p w 3.4 (the corresponding chemical potential 
fi —0.2) of the particles the antisymmetric state be- 
comes stable; see the spectrum of the stable state shown 
in Fig.[5][d). The peculiarity of this state stems from the 
fact that the ground state is strongly depleted in the areas 
with repulsive nonlinear interatomic interaction. There- 
fore, for large negative chemical potentials some excita- 
tions can be considered as motion of weakly interacting 
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drops of the condensate and this interaction becomes ex- 
ponentially weak for large negative fj,. It means that the 
eigenvalues appear in the spectrum with very small ab- 
solute values and the chain of the droplets becomes very 
flat so even relatively weak external action results in the 
significant motion of the mutual phases of the droplets. 
Therefore, in this case even relatively weak external per- 
turbation can result in the strongly nonlinear dynamics 
and in the decay of the condensate. As a result, even rel- 
atively weak perturbations can destroy stability due to 
nonlinear effects, while the state remains linearly stable. 

We now consider how the stability of the antisym- 
metric state changes when vq deviates from zero and 
the ground state has nonzero current. In general, the 
period of the coefficients is equal to 2tt; therefore, the 
width of the zone is equal to 1. However, as mentioned 
above, at vq = the period of the coefficients in the lin- 
earized equation is actually equal to ir; - see FigJHJa). 
For nonzero vq the period becomes equal to 2ir and gaps 
appear in the spectrum, as one can see from Figs. HJb) 
and[B][c). One can also see that the symmetry k — > — fc 
becomes broken for nonzero current. For relatively small 
values of vq the background remains stable and thus one 
can claim existence of stable ground states with nonzero 
current. 

At some intermediate value of vq the modes start to 
merge, generating instability. The onset of the instabil- 
ity goes to the Josephson critical current with the deple- 
tion, but for finite depletion the instability appears at 
the current lower the critical one. The development of 
the instability is illustrated in Fig. SJ We took the ini- 
tial condition in the form of the numerically found sta- 
tionary solution perturbed by weak noise and performed 
the modeling. One can see that the noise destroys the 
state during relatively short time. At vq = 1 the states 
transform into the symmetric mode considered above. As 
noted, this state is unstable. The gaps at the boundary 
of the zone become closed because the period of the co- 
efficients becomes equal to it. 

Finally, we remark that the fact that the obtained spec- 
tra are of Bogoliubov type means that a perturbation 
moving (with subcritical velocity) with respect to the 
condensate does produce scattering (a certain amount of 
scattering may be seen, however, shortly after the intro- 
duction of the perturbation). This applies to the case of a 
perturbation that is resting with respect to the nonlinear 
potential, a problem to be addressed in the next section. 
The more general case of perturbations moving with re- 
spect to both the nonlinear lattice and the condensate 
will be considered elsewhere. 



IV. SUPERFLUIDITY OF BOSE-EINSTEIN 
CONDENSATE 

It is well known that the scattering of the particles 
can be interpreted as Cherenkov synchronism between an 
obstacle and the propagating eigenmodes of the electron 




(b) k (d) 
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FIG. 6: (Color online) Sound wave spectra for the state with 
H = — 2 and the quasivelocities vo = (a), vo = 0.33 (b), 
vq = 0.66 (c) and vo = 1 (d). The antisymmetric state (a) 
and the symmetric one (d) modes correspond curves 2 and 1 
in Fig. [1] respectively. Blue curves correspond to imaginary 
part u) of the eigenvalue (left vertical axis) and the red curves 
to the real part 7 of the eigenvalues (right vertical axis). 



gas. When the system is described by the Schrodinger 
equation, then a localized obstacle can be represented as 
localized linear potential in the equation. For the lin- 
ear Schrodinger equation, with the dispersion of the el- 
ementary excitations parabolic, the synchronism with a 
moving obstacle is inevitable unless the obstacle is rest- 
ing with respect to the condensate. However, in the case 
of nonlinear Schrodinger equation with an uniform re- 
pulsive (defocusing) nonlinearity the dispersion becomes 
linear at low k, implying the existence of a critical ve- 
locity. The scattering occurs only if the velocity of the 
obstacle is larger than the critical velocity. 

In the previous section we have shown that there exist 
the spectra of linear excitations in Bose- Einstein conden- 
sates with spatially inhomogeneous nonlinear interaction 
of Bogoliubov kind with nonzero v g at k — 0. There- 
fore, one can expect that at least an obstacle resting rel- 
atively to the nonlinear periodical potential will not lead 
the scattering of the Bose- Einstein condensate flowing in 
the system (the condensate is flowing relatively to both 
the periodical potential and the obstacle). To confirm 
this hypothesis we performed numerical simulations of 
the condensate flowing in the large annular system with 
periodic nonlinear potential. The results are summarized 

Fig. m 

One can see that for the stationary state the current is 
evenly distributed in space. If an obstacle leads to scat- 
tering, then the distribution of the current becomes in- 
homogeneous and, after some time, the averaged current 
tends to zero. We always calculated the distributions 
of currents assuming the initial condition the antisym- 
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metric ground state discussed above and perturbing the 
equation by localized linear potential resting in respect 
to the periodical nonlinear potential. We carried out the 
simulations for very long time and no continuous scat- 
tering has been observed; the dependence of the current 
averaged over space is shown in Fig.[71[c). 

The results of the numerical simulations indicate that 
introduction of the obstacle leads to very small perturba- 
tion of the current and the perturbation does not grow 
in time. This phenomenon is not sensitive to the rela- 
tive position of the obstacle in respect to the nonlinear 
periodical potential; see Fig. [7] We also verified that the 
perturbation of the current does not grow at arbitrary 
point of the system. For comparison, we show in the 
bottom part of Fig. [TJc) the dependence of the averaged 
current in the system with the same average particle den- 
sity and the same initial current but without nonlinear 
interatom interaction. It is evident that the current oscil- 
lates and that the calculation time is considerably longer 
than the characteristic scattering time. The oscillation 
of the current in the linear system is, of course, periodic. 
Therefore, we can conclude that superfluidity of the BEC 
is observed in the system with periodic nonlinear poten- 
tial. 



V. CONCLUSION 

We have explored the combined effect of nonlinearity 
and periodicity of nonlinear inter-atomic interaction in 
the toroidal traps where the scattering length is periodi- 
cally modulated along the ring. First, we systematically 
studied periodic matter waves for both stationary states 
and states carrying nonzero current. The solutions found 
numerically have been parametrized by a chemical po- 
tential. We have shown for stationary states that when 
the nonlinear interaction changes its sign the condensate 
exists in the form of interacting drops localized in the 
areas with attractive interaction which rather counter- 
intuitively implies the existence of the condensate with 
negative chemical potential despite of the dominating po- 
tential being repulsive. We also studied numerically the 
states with non-zero current and their relation to the sta- 
tionary states. We have shown that the current depends 
periodically on the phase gradient and displays a strong 
fragmentation for large negative values of the chemical 
potential when the condensate is significantly depleted 
in the areas with repulsive interatomic interaction. We 
interpreted the fragmented condensates as Josephson n- 
j unctions and the presence of the persistent current as a 
BEC analog of dc- Josephson effect which in contrast to 
a standard Josephson junction occurs in the system with 
and attractive effective inter-atomic potential. 

Simultaneously, we studied sound waves propagating 
against the background. By inspecting the spectra of 
linear excitations on the ground states we investigated 
the stability of the condensate and we have shown that 
the spectrum are of Bogoliubov type. We found that the 
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FIG. 7: (Color online) The dynamics of BEC in the 
presence of an obstacle described by the potential V = 

( — 

clq exp(— (x x °' ). Panels (a) and (b) show the relative po- 

sitions of the nonlinear periodical potential and the linear 
potential modeling an obstacle, xo = for panel (a) and 
xo — tt/2 for panel (b). The dynamics of the averaged current 
(J) is shown in the upper part of panel (c) for the relative po- 
sition of the potentials shown in panels (a) (black line) and 
(b) (red line). The dynamics of the averaged current for the 
case without nonlinear interatom interaction is shown in the 
bottom part of panel (c). The averaged density and the ini- 
tial velocity of the particles are the same as for the case with 
nonlinear interactions. Parameters used are jj, = —2, vq = g, 
a = 0.15, wo = 0.8, and L = 64/7T. Within the time scale 
shown in the figure the condensate makes about three rounds 
in the trap with a duration of each round of 512-7T in dimen- 
sionless time units. 



ground states consisting of interacting drops can be sta- 
ble, which provides the strong evidence that these states 
and dc- Josephson effect can be experimentally observed. 
On the basis of numerical simulations we have demon- 
strated that superfluidity of Bose-Einstein condensate 
may take place in the system with inhomogeneous nonlin- 
ear interatom interaction. In particular, we have shown 
that a localized defect does not lead to the scattering of 
the condensate when the defect is at rest with respect to 
the nonlinear potential while the condensate is moving 
in respect to both the periodic potential and the defect. 
We have observed this behavior in the systems with both 
an effective repulsive and attractive interatom interac- 
tion. The case when the defect is moving in respect to 
both the lattice and the condensate will be considered in 
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a separate paper. 
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